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ABSTRACT 



In this paper we compute the emission coming from the direct recombination of free electrons to a given shell (n > 2) during the 
epoch of cosmological hydrogen recombination. This contribution leads to a total of one photon per recombined hydrogen atom and 
therefore a ~ 30 - 88 % increase of the recombination spectrum within the frequency range 1 GHz < v < 100 GHz. In particular the 
Balmer-continuum emission increases the distortion at v ~ 690 GHz by ~ 92%. With our 100 shell calculations for the hydrogen atom 
we find that a total of ~ 5 photons per hydrogen atom are emitted when including all the bound-bound transitions, the 2s two-photon 
decay channel and the optically thin free-bound transitions. Since the direct recombination continuum at high n is very broad only a 
few n-series continuua are distinguishable and most of this additional emission below v < 30 GHz is completely featureless. 
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1. Introduction 

With the advent of accurate observations of the Cosmic 
Microwave Background (CMB) temperature and polariza- 
tion anisotropies it becomes increasingly important to un- 
derstand the dynamics of cosmological recombination on 
the percent level precision . Seve ral authors tLeung et alJ 
2004; Dubrovich & Grachev I20051 IChluba & Sunvaevl 



Kholupenko & Ivanchik 2006) have discussed physical pro- 
cesses leading to perc ent level correct ions to the results of the 
standard computation ( Seage r et alJ200 0) for the ionization his- 
tory. Recently, it has been shown that also the treatment of 
the populations in the angular momentum sub-states of hydro- 
gen has percent level impact on the ionizatio n history, espe- 
cially at redshifts z < 800 - 1000 dRubino-Martin et alJl2006l 
IChluba etal1l2006L hereafter RMCS06 and CRMS06 respec- 
tively) and it is expected that more physical processes altering 
the dynamics of c osmological recombination may be realized. 
iLewis etaf] J2006) have made some first steps towards quanti- 
fying the possible impact of percent level corrections to the ion- 
ization history on the estimation of cosmological parameters. 

In the future it may become possible to directly observe 
the spectral distortions of the CMB arising during the epoch 
of recombination. This would in principle open an alterna- 
tive way to determine cosmological parameters like the baryon 
and total matter density. Several authors have discussed the 
distortions arising due to the hydrogen higher level b ound 
bound transitio ns ( DubroTjcl^l9J75yLhibarskii &^unvaevll983 



Rvbicki & dell' Antonio! 119931: iD ubrovich & Stolvarov 1995 



Dubrovich & Shakhvorostova 2004; Kholupenko et al. 2005 
Wong et alJ2006L RMCS06~CRMS06'). 

In the standard computations within the cosmological con- 
text the direct recombination to the ground state of hydrogen 
is neglected, since this transition is so optically thick during 
the whole epoch of recombination that escape of photons in 
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the Lyman-continuu m is considered impossible (Zeld ovich et alJ 
1968: IPeebiesll968l) . On the other hand electrons can reach any 
of the other levels (with n > 2) by direct recombination and the 
emitted continuum photons escape freely, just because after their 
creation they will not encounter another optically thick transition 
within the hydrogen atom. 

Every successful recombination (i.e. when the electron is 
reaching the ground-state and remains there) will therefore lead 
to the emission of at least two photons, one from the direct re- 
combination and (because direct recombinations to the ground- 
state are neglected) at least one within the cascade to the ground- 
state. Here we show that for hydrogen a total of ~ 4 photons per 
neutral hydrogen atom is produced within the bound-bound cas- 
cade and 2s two-photon decay channel. Therefore, one expects 
a ~ 25% addition to the total number of emitted photons during 
the epoch of recombination due to the direct recombination to 
states with n > 2. 

In this paper we discuss the hydrogen recombination spec- 
trum in the frequency range 100 MHz < v < 3000 GHz, with 
special emphasis on the contribution from the direct recombina- 
tion lines (for more details on the bound-bound emission within 
this context see RMCS06 and CRMS06). We use the solution for 
the recombination history and evolution of the hydrogen popula- 
tions as obtained in CRMS06. In those computations a maximum 
of 100 shells following the populations of all the angular mo- 
mentum sub-states separately and including /- and n-changing 
collisions was treated. We refer the reader to this paper for more 
details on the computations. 



2. Spectral distortion due to direct recombination 

In order to compute the spectral distortion arising due to direct 
recombinations to a given level (n, /) one has to consider the 
emission and absorption of photons due to this process. The the 
effective rate for the change of the population, TV;, of a hydrogen 
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level i due to recombinations from the continuum is given by 



dt 



= N e N p R ci - NiR ic , 



(1) 



where N e and N p are the free electron and proton number den- 
sities, R ci is the corresponding recombination rate and R ic the 
photoionizationrate. Using the definition of R C j and R, c in terms 
of the photoionization cross section, cr IC (v), one can write the ef- 
fective change of the photon number density, N v = I v /hv, due to 
direct recombination to level i at frequency v as 



1 8Ny 

c ~dt 



N e N p MT e )<r k 



2v 2 



N iC r ic N v . (2) 



Here stimulated recombination is included. From the Saha- 

(\LTE 3/2 
m) = Kot) eEllkTc > 

where Ej is the ionization energy, T e is the electron tempera- 
ture, which always is very close to the radiation temperature 
T y = 7o(l + z) with T = 2.725 K, and gi = 2(2/ + 1) is the 
statistical weight of the level. It is clear from equations Q and 
(0 that the total change of the number of photons has to be iden- 
tical to the number of recombined electrons. 

Now, given the solution for the recombination history and 
evolution of the populations Ni, one can obtain the solution for 
the change of the radiation field in the optically thin limit as- 
suming a pure blackbody ambient photon field. We find for the 
spectral distortion of the CMB at observing frequency v due to 
direct recombination to level i at redshift z — 



-s,r 



'cNiQ- ic (v z ) 
H(z)(l+z) 



N e N „ 



■MT e )e* 



-m _ i 



dz, 



(3) 

with v, = v(l +z) and where H(z) is the Hubble-expansion factor, 
1 + Zt = Vfc/v corresponds to the redshift at which the emission 
and photoionization threshold frequency, v !C , are equal and B v is 
the CMB blackbody spectrum today. 

In the limit of n » 1 the ionization threshold corresponds to 
hv„ = x/n 2 an d the energy of the n a transition is hv„ a = 2^/n 3 , 
i.e. n/2 times smaller. Therefore in our computations with « max = 
100 at z — the lowest frequency we reach is expected to be of 
the order vi ow ,ft> ~ 200 - 300 MHz instead of vi ow , bb -4-6 MHz 
for the bound-bound transitions . The width of the recombina- 



tional line is of the order /iAv„ ~ kT^ k , where T t 



^peak 



is the temperature of the electrons at the redshift, where the main 
contribution to the direct recombination emission for shell n ap- 
peared (typically zJJ eak ~ 1300). Hence, one expects that for 
hv„ <sc kTl^ the contributions to the free-bound continuum will 
become very broad. For the Balmer and Paschen continuua one 
has hv2 - lOkT^^ and /1V3 ~ 5 kT^ k , so that the contribution 
due to these transitions will be narrow, while for n » 1 they will 
be very broad (see Fig.0. 



3. Results and discussion 

In Fig. [2 we present the full hydrogen recombination spec- 
trum including the emission due to bound-bound transitions, the 
2s two-photon decay and direct recombinations to shells with 
n > 2. At high frequencies the free-bound emission shows nar- 
row features corresponding to the Balmer, Paschen, Brackett and 
Pfund continuua. At low frequencies the free-bound continuua 
overlap very strongly and the total continuum emission becomes 
completely featureless. The inlay plot in the upper panel shows 



the free-bound continuum emission for three well separated high 
levels. These contributions are very broad and each of them lies 
by a factor of ~ 10 below the total free-bound emission spec- 
trum. The total level of the emission is only reached after sum- 
ming over many shells. As the inlay plot in the lower panel il- 
lustrates the sum of all contributions (here mainly bound-bound 
and free-bound) still has a significant modulation ranging from 
~ 4-35% in the frequency band 1 GHz < v < 30 GHz. The slope 
of the free-bound distortion for 1 GHz < v < 10 GHz is ~ 0.6 
and therefore slightly steeper than for the contribution from the 
bound-bound transitions alone (~ 0.46). At v < 1 GHz one still 
expects an increase of emission when including more than 100 
shells. 

In Fig.^we also show the free-bound contribution for differ- 
ent values of n max < 100. From this one can conclude that within 
our assumptions the free-bound emission spectrum is converged 
to a level of better than 1% at v > 5 - 10 GHz and better than 
10% for 1 GHz < v < 2 GHz. We have also checked how the 
free-bound emission is depending on the treatment of the an- 
gular momentum sub-states and found that at high frequencies 
the solution, when assuming full statistical equilibrium within 
the shells for n > 2, is not very different. At low frequencies, 
like in the case for bound-bound transitions (see CRMS06), the 
difference is much larger, again showing that it is important to 
follow the populations of all the angular momentum sub-states 
separately. 

Figure [2] shows the total number of photons per hydrogen 
atom emitted during recombination for a given n-series (sum of 
all bound-bound transitions to one shell with fixed n) and its con- 
tinuum (the corresponding bound-free contribution). In the con- 
sidered case, one can see that for n > n CI - 30 the total emission 
is dominated by the contribution from direct recombinations. For 
n < n CI the contribution from bound-bound-transitions domi- 
nates. This behavior shows that for the lower shells cascading 
electrons are very important, whereas for larger n the electrons 
reach a given shell mainly via direct recombinations. Obviously 
including more than 100 shells will lead to an increase of n cr , but 
one does not expect a significant addition to the total number of 
emitted photons. There is a strong difference in both curves when 
the computation is done assuming full statistical equilibrium for 
shell with n > 2. In this case the curves intersect at « cr ~ 59 
instead of n cl ~ 30. This shows that following all the angular 
momentum sub-states the lack of strong redistribution to states 
with / » 1 disfavours the cascade with transitions An <k n. 

Knowing the asymptotic behavior of the distortion at low 
frequencies, AI V = A [v/lGHz]^ lO^JirT 2 s" 1 Hz" 1 sr 1 , from 
Fig [2 one can find the number of photons emitted per hydro- 
gen atom (present density A^h = 1.9 x 10~ 7 cirT 3 ) at v < vq, 
i.e. N y /N H = ^ J V0 % dv = 3.3 x 10~ 2 f [v /lGHzf. Within 
the range 1 GHz < v < 10 GHz we obtained Ajf -2.2 and 
/?bb ~ 0.46 for the bound-bound and A* ~ 0.59 and /fo ~ 0.6 for 
the free-bound distortion. Considering the bound-bound distor- 
tion assuming that the <?-transitions give the main contribution 
and that most of the photons are released at z e m ~ 1300, with 
hv^ = 2 X /z em n\ one finds Nf{n > n )/N H = 3.3 x 10- 2 [6.6 x 

10 6 /Zem]' 3 j « 3/i ~ 8.0 n 138 and | dA^ b / dn,,l ~ 1 1 « 238 A^h for 
the contribution of the whole series for given shell «o- Similarly, 
for the free-bound emission with /iv° bs = ;f/z e m« 2 one has 
A^(n > no)/N H ~ 3.6 n L2 and |dA#7dn | ~ 4.3 n 2 2 N H , al- 
though for no » 1 due to the large width of the free-bound con- 
tributions one expects this estimate to be rather crude. We found 
that IdA^/dnol ~ 2.5n L9 A^e represents the result shown in 
Fig.|5]for large no very well. One should mention that for n » 1 
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v [GHz] 

Fig. 1. The full hydrogen recombination spectrum including the free-bound emission. The results of the computation for 100 shells 
as presented in CRMS06 were used. The contribution due to the 2s two-photon decay is also accounted for. The dashed lines indicate 
the expected level of emission when including more shells. In the upper panel we also show the free-bound continuum spectrum 
for different values of n max (dashed-dotted). The inlay gives the free-bound emission for n = 30, 40 and 50. The lower panel shows 
the distortion relative to the CMB blackbody spectrum and the inlay illustrates the modulation of the total emission spectrum for 
1 GHz < v < 30 GHz in convenient coordinates. 



4 



Chluba and Sunyaev: Free-bound emission from cosmological hydrogen recombination 



IFI 1 ' 1 1 1 ' 1 ' 1 <-i 




n (principal quantum number) 



Fig. 2. Total number of photons per hydrogen atom released dur- 
ing recombination for a given n-series and its continuum for 
"max = 100. Also the results assuming full statistical equilibrium 
within the shells for n > 2 (« sp ii t = 2) is presented. 

Table 1. Total number of photons per hydrogen atom emitted 
during recombination within a given series and its continuum. 
Adding all the contributions one obtains 4.98 photons per hy- 
drogen atom. 



Series 


n 


bound-bound 


continuum 


sum 


Lyman 


1 


4.28 x 10" 1 





4.28 x 10 1 


B aimer 


2 


7.83 x 10-' 


2.11 x 10-' 


9.94 x 10~' 


Paschen 


3 


4.91 x 10-' 


1.37 x 10-' 


6.28 x 10-' 


Brackett 


4 


2.38 x 10-' 


9.59 x 10~ 2 


3.34 x 10-' 


Pfund 


5 


1.62 x 10-' 


7.22 x 10~ 2 


2.34 x 10-' 


All 


1-100 


2.84 


1.00 


3.84 


2s-decay 


1 





1.14 


1.14 



signal under discussion is close to 2 x 10~ 7 B V and still has vari- 
ability with well defined frequency dependence on the level of 
several percent (Fig.[0. Also in this band (A < 21 cm) one does 
not expect other sources with similar frequency dependence. 

Acknowledgements. We wish to thank J. A. Rubino-Martin for useful discussion 
on this problem at the initial stage. 
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and | An | <k n induced recombinations and stimulated transitions 
play an important role. 

In Tabled we give the values of the total number of photons 
per hydrogen atom emitted in the B aimer, Paschen, Brackett and 
Pfund-series and their continuua. Also we listed the total num- 
ber of emitted photons in all n-series and their continuua and the 
contribution from the 2s two-photon decay. As expected the total 
number of emitted photons due to free-bound transitions (within 
the accuracy of our calculation) is identical to the total number 
of hydrogen atoms. The number of photons emitted within the 
Balmer-series and B aimer-continuum is also very close to one 
photon per hydrogen atom. One can see that roughly 2.8 pho- 
tons per hydrogen atom are released within the bound-bound 
transitions. For the free-bound distortion 80% of the photon are 
emitted for n < 14 - 15, 90% for n < 26 and 95% for n < 41. 
On the other hand within the bound-bound transitions one finds 
80% for n < 6 - 7, 90% for n < 1 1 - 12 and 95% for n < 18. 
Focusing on the Lyman-series and the 2s two-photon decay con- 
tribution one can see that ~ 43 % of all electron go through the 
Lyman-a and ~ 57 % of electron reach the ground state via the 
2s two-photon decay channel. During the epoch of hydrogen re- 
combination a total of ~ 5 photons per hydrogen atom are pro- 
duced. Therefore, recombination of hydrogen slightly increases 
the specific entropy of the Universe (photons perbaryon). 

Observations of the recombination features at frequencies 
v > 1412 MHz might become feasible since the strength of the 



